function q = Elnqmu(mix)

% E[ln q(mu)]

[K D] = size(mix.centres);
Dlog2pi = D*log(2*pi);

T = mix.varposterior.T;

q = 0;
for k = 1:K
    q = q - 0.5*D - 0.5*Dlog2pi - 0.5*logdet_chol(T(:, :, k));
end
